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^sO | We propose a fast algorithm for computing the expected tranche 

' loss in the Gaussian factor model with arbitrary accuracy using Her- 

■ mite expansions. No assumptions about homogeneity of the portfolio 

! are made. It is intended as an alternative to the much slower Fourier 

' transform based methods 

a" 

> 1 The Gaussian Factor Model 

X. 

Let us consider a portfolio of N loans. Let the notional of loan i be equal to 
the fraction of the notional of the whole portfolio. This means that if loan 
i defaults and the entire notional of the loan is lost the portfolio loses fraction 
fi or 100/1% of its value. In practice when a loan i defaults a fraction r« of 
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its notional will be recovered by the creditors. Thus the actual loss given 
default (LGD) of loan % is 



LGDi = fi(l — n) (1) 

fraction or 

LGDi — 100/j(l - rj)% (2) 

of the notional of the entire portfolio. 

We now describe the Gaussian m-factor model of portfolio losses from 
default. The model requires a number of input parameters. For each loan i we 
are give a probability pi of its default. Also for each i and each k — 1, . . . , m 
we are given a number w it k such that J2t=i w ik < 1- The number w it k is 
the loading factor of the loan i with respect to factor k. Let (pi, ... , (f> m and 
4>\i = 1, . . . , N be independent standard normal random variables. Let 
be the cdf of the standard normal distribution. In our model loan i defaults 
if 



k=l 



\ 



k=l 



This indeed happens with probability p^. The factors cf>i, . . . , (j) m are usually 
interpreted as the state of the global economy, the state of the regional econ- 
omy, the state of a particular industry and so on. Thus they are the factors 
that affect the default behavior of all or at least a large group of loans in the 
portfolio. The factors 1 , . . . , <p N are interpreted as the idiosyncratic risks of 
the loans in the portfolio. 
Let Ii be defined by 

li I{loan i defaulted} (4) 

We define the random loss caused by the default of loan i as 

U = fi(i - n)l h (5) 

where r; L is the recovery rate of loan i. The total loss of the portfolio is 

£ = E^i (6) 

i 

An important property of the Gaussian factor model is that the Lj's are 
not independent of each other. Their mutual dependence is induced by the 
dependence of each Lj on the common factors <pi, . . . , <p m . Historical data 
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supports the conclusion that losses due to defaults on different loans are 
correlated with each other. Historical data can also be used to calibrate 
the loadings u^.the Lj's are not independent of each other. Their mutual 
dependence is induced by the dependence of each Lj on the common factors 
0i, ... , <f> m . Historical data supports the conclusion that losses due to defaults 
on different loans are correlated with each other. Historical data can also be 
used to calibrate the loadings w^. 



2 Analytic Approximation to the Joint Dis- 
tribution of </>!,..., (j) m and L 

When the values of the factors </>i, . . . , m are fixed, the probability of the 
default of loan i becomes 

i = 3-1 ( Pi- E k w^ k \ 

The random losses Lj become conditionally independent Bernoulli vari- 
ables with the mean given by 

E cond (L i ) = f i (l-r i )p i (8) 

and the variance given by 

VAR cond (Li) = if (1 - n) V(l - V 1 ) (9) 

By the Central Limit Theorem the conditional distribution of the portfolio 
loss L given the values of the factors (pi, . . . , <p m can be approximated by the 
normal distribution with the mean 

Econd(L) = y~] Econ^Lj) (10) 

i 

and the variance 

VAR cond (L) = VAR cond {U) (11) 

i 

Then the joint distribution of the factors 0i, . . . , <p m and the portfolio loss 
L can be approximated by a distribution with density 

m 

p(01, . . . ,(f) m ,L) = p G ,E cond (L),VAR cond (L)( L ) II PG,0,l{<Pk), (12) 

fc=l 
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where PG,m,v(%) stands for the Gaussian density with mean m and variance 
v. 



3 Expected Loss of a Tranche of Loan Port- 
folio 

Let < a < b < 1. We define a tranche loss profile Tl a ^(x) by 

_ min(b - a, max(x - a, 0)) 
J- i a ,b{x) — 7 U-jJ 

o — a 

Number a is called the attachment point of a tranche, while b is called the 
detachment point of a tranche. The expected loss of a tranche is then 

TLoss(a, b) = J Tl a , b (L)p(<p u ...Am, L)d$ x . . . <p m L (14) 

This can be rewritten as a double integral 

„ „ TO 

TLoss(a, b) = J J Tl a ^ b {L)p GtEcond{L) y ARcond{L) {L)dL J] p G ,o,i(0fc)^i • • • <Pm 



k=i 



(15) 

The inside integral with respect to L is easily done analytically because L 
has a simple normal distribution for fixed values of the factors <pi, . . . , (j> m . 
The outside integral has to be computed numerically. However, since it is 
an integral of a bounded smooth function with respect to m-dimensional 
Gaussian density, it is one of the simpler integrals to compute numerically. 



4 Numerical Example 

In this section we apply the proposed algorithm to the single factor Gaussian 
model of a portfolio with 125 names. We choose a 125 name portfolio because 
it is the size of the standard DJCDX.NA.IG portfolio. We choose a single 
factor model because it is the one most frequently used in practice. We 
evaluate the expected loss for four different tranches. All tranches have 
attachment point a = or a = 0%. The detachment points are 3%, 7%, 10% 
and 15%. We take these detachment points because they are the ones most 
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frequently used in practice in order to evaluate the base correlation. The 
parameters of the porfolio are 



h = ^ 



1 

125 



Pi 



0.015 + 
0.5- 



0.05(z 



124 

o.Kz-1) 



Wil 



0.5 



124 

o.iU- v 



124 



(16) 



In Figured] we compare the expected loss computed using 10 6 Monte Carlo 
samples with the expected loss computed using formula (JT5j) . The agreement 
between the two is good. 

Figure 1: Loan Portfolio Tranche Loss in the Gaussian Single Factor Model 
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To obtain the results in Figure ^ we only needed to perform a single one 
dimensional numerical integration for each tranche. This is an improvement 
over the Moody's FT method j2] which requires computing a large number 
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of Fourier transforms for each tranche. Each individual Fourier transform is 
as computationally expensive as (fT5j) . 



5 The Hermite Expansion of the Conditional 
Distribution of the Portfolio Loss L 

When the approximation to the conditional distribution of L given by the 
Central Limit Theorem is deemed insufficiently accurate, an arbitrarily ac- 
curate representation of the conditional distribution of the portfolio loss L 
can be obtained from its Hermite series expansion. For historical reasons this 
expansion is also known as the Charlier series expansion [S], 0. 

Let F(x) be the c.d.f. of the conditional distribution of the portfolio loss 
L. So that 

P(L <x) = F(x) (17) 

For each fixed value of the factors (pi, ... , <p m we define the normalized con- 
ditional loss L by 

L — L ~ E «U L ) (18) 
/VAR cond (L) 



Let F(x) be the c.d.f. of the distribution of the normalized conditional 
portfolio loss L. So that 

P(L < x) = F(x) (19) 
We define the Hermite polynomial H n (x) of degree n by 

H n (x) = (-ire^-^e^ (20) 

Let c n be defined by 



ni J-oo 

Then we have 



H n (x)dF(x) (21) 



Hx) = E r *Hi{t)ZjLdt (22) 
i=o J -°° v2tt 
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The series above converges in the sense of distributions (generalized func- 
tions) 0] . A good reference on the theory of distributions (generalized func- 
tions) is [3] . Let us pick a finite N. Then we have 

f H i(t)^=dt (23) 
i=o J -°° v2tt 

As before the approximation is in the sense of generalized functions. Equation 
(I23|) implies that the distribution of the normalized conditional portfolio loss 
L can be approximated by a distribution with the density 

~p( x ) = J2c l H t (x)^-L (24) 

,-_n V27T 



The function p(x) is not necessarily nonnegative and therefore may not be 
a probability density in the strict sense. However, as is explained in jl], this 
does not affect the validity of our final result (jl5j) . Therefore we may treat 
p(x) as a real probability density. 

The distribution of the unnormalized loss L can be approximated by a 
distribution with density 

/ x - E cond( L ) 



P ( X ) = V y Ci H< " - Econd(L) 6 - ^ (25) 



o ^VAR cond (L) \ AR cond (L) ) V2ir 

The joint distribution of the factors 0i, . . . , <p m and the portfolio loss L 
can be approximated by a distribution with density 

m 

Pjoint{4>l, ...,(j)m,L) = p(L) Yl PG,0,l{<Pk), (26) 

fe=l 

where Pg,o,i( x ) stands for the Gaussian density with mean and variance 1. 

Observe that the coefficient c n depends only on the moments of the dis- 
tribution F(x). Since Lj's are independent Bernoulli random variables these 
moments are known analytically. Thus in the case under consideration all 
the known analytically. 

If in equation (|23j) we set iV = 1 we obtain the normal density approx- 
imation proposed in the previous section. We show later that it gives good 
numerical results even when the portfolio size is too small for the normal 
approximation to be accurate. 
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6 Expected Loss of a Tranche of Loan Port- 
folio 



The expected loss of a tranche can be written as a double integral 

TLoss{a, b) = J I Tl ajb (L)p(L)dL J[ PG,o,i(0fc)#i. ■■■4>m (27) 

k=l 

The inside integral with respect to L can be done analytically for fixed values 
of the factors 0i, . . . , <fi m . The outside integral has to be computed numer- 
ically. However, since it is an integral of a bounded smooth function with 
respect to m-dimensional Gaussian density, it is one of the simpler integrals 
to compute numerically. 



7 Numerical Example 2 

In this section we test the proposed algorithm on several portfolios of smaller 
size. For these portfolios the approximation to the conditional distribution of 
the portfolio loss L given by the Central Limit Theorem is not very accurate, 
because of their small size. However, the Hermite expansion produces very 
good results. We apply the proposed algorithm to the single factor Gaussian 
model of a portfolio with n names. We take n to be 25 (size of DJ iTraxx 
Australia), 30 (size of DJ iTraxx ex Japan), 50 (size of DJ iTraxx CJ) and 
100 (size of DJCDX.NA.HY). We choose a single factor model because it is 
the one most frequently used in practice. For each n we compute the loss 
of the equity tranche with the attachment point a = or a = 0% and the 
detachment point 3%. The parameters of the porfolio are 

/, = i 

n 

p, = 0.015 +O-O50-D 



n — 1 

0.1U - 1) 

n = 0.5 K — -L 

n — 1 

0.1(2-1) / \ 

wn = 0.5 y — A 28 

71—1 



where % = 1, . . . , n. Finally, we choose iV = 5 in f!23l) . 
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In Figure 121 we compare the expected loss computed using 10 6 Monte 
Carlo samples with the expected loss computed using formula §I7\ . 1 The 
agreement between the two is good. 



Figure 2: Equity Tranche Loss in the Gaussian Single Factor Model 
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To obtain the results in Figure |2l we only needed to perform a single one 
dimensional numerical integration for each tranche. This is an improvement 
over the Fourier transform based methods j2] which require computing a 
large number of Fourier transforms for each tranche. Each individual Fourier 
transform is as computationally expensive as (|27j). 

The Hermite expansion ()23j) can be used to achieve arbitrary accuracy 
when the normal approximation is insufficiently accurate. The proposed 
algorithm is fast because the inside integral in (|2*7|) can be done analytically. 

We also comment that the algorithm can be extended trivially to the case 
of non-constant recovery rates and recovery rates correlated with the state 
of the factor variables. 

1 The author has the code implementing the algorithm described here in MATLAB, 
VBA for Excel and C. 
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